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A  NEW  METHOD  OF  COMPUTING  PENETRATION  VARIABLES  FOR 
SHAPED-CHARGE  JETS 


In  earlier  attempts  to  describe  the  penetration  of  shaped-charge  jets, 
the  total  depth  of  penetration  was  computed  from  J"*  Udt.  Evaluation  of  the 
integral  has  not  been  possible  because  the  penetration  velocity  U  for  the 
discontinuous  portion  of  the  jet  has  not  been  properly  defined.  In  the 

following  analysis,  it  is  postulated  that  the  total  penetration  can  be 

c  \ — * 

represented  by  j  Udt  +  y  AP,  where  the  integral  is  taken  over  the  penetration 
produced  by  the  continuous  portion  of  the  jet  and  the  summation  is  taken  over 
the  penetration  produced  by  the  discrete  particles  resulting  from  jet  breakup. 


It  is  shown  that  numerical  values  of  the  penetration  variables  for  the  lOpmm 
unconfined  test  charge  are  compatible  with  the  modified  penetration  equation. 
Scaling  relations  are  introduced,  and  the  penetration  equations  are  given  in 
size  Independent  variables,  which  are  substantiated  experimentally.  Although 
it  was  assumed  that  the  incompressible  approximation  provides  an  adequate 
relationship  between  the  jet  and  penetration  velocities  for  a  continuous  jet. 
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materials  do  not  seriously  alter  this  relation.  The  incompressible  approximation 
does  underestimate  the  pressure  by  an  appreciable  amount. 
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INTRODUCTION 


A  shaped  charge  is  a  jet-forming  explosive  device  containing  an  axially 

symmetric  cavity  lined  with  a  thin  metal  liner.  The  geometrical  design  of  a 

shaped  charge  can  be  varied  over  a  wide  range;  however,  cylindrical  charges 

containing  conical  cavities  have  been  most  extensively  investigated.  Line 

drawings  illustrating  the  jet  formation  and  penetration  process  for  a  typical 

shaped  charge  are  presented  in  Figure  1.  When  the  explosive  is  detonated,  the 

high  pressure  of  the  detonation  products  collapses  the  metal  liner  upon  itself 

forming  a  high-velocity,  forward-moving,  metallic  jet,  which  is  able  to 

penetrate  high-strength  materials  such  as  armor.'1'  Different  elements  of  the 

liner  experience  different  collapse  velocities  and  the  resulting  jet  contains 

a  spectrum  of  velocities  with  forward  portions  of  the  jet  traveling  faster 

2  3 

than  succeeding  portions.  ’  The  velocity  gradient  in  the  jet  causes  it  to 
lengthen  as  it  travels  forward  from  the  charge.  Initially  the  jet  lengthens 
in  a  ductile  manner,  but  eventually  the  ductility  of  the  material  is  exceeded 
and  the  jet  fractures  into  a  series  of  discrete  particles. 

The  pressure  exerted  by  the  impinging  jet  greatly  exceeds  the  yield 
strength  of  the  target  material,  and  penetration  takes  place  by  a  hydrodynamic 
process  as  illustrated  in  Figure  1.  Theories  describing  the  penetration  of  a 

shaped- charge  jet  were  advanced  independently  by  Pugh'1'  in  the  United  States  and 

* 

by  Hill,  Mott,  and  Pack  in  Great  Britain.  Both  theories  were  based  on  steady- 
state  incompressible  fluid  flow  and  resulted  in  an  equation  of  the  form 

\Pj(Vj  -  U)2  =  p^U2,  (1) 

where  p.  is  the  jet  density;  p, ,  the  target  density;  V.,  the  velocity  of  the 
J  j 

jet  at  the  stagnation  point;  U  =  dp/dt,  the  penetration  velocity;  and  a 

parameter  introduced  to  account  for  differences  in  the  transfer  of  momentum 

to  the  target  from  a  continuous  jet  and  from  a  jet  of  discrete  particles.  The 

accuracy  of  the  penetration  equation  has  been  extensively  investigated  by 


The  existence  of  independent  work  by  Hill,  Mott,  and  Pack  was  cited  in 
reference 
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Eichelberger  ,  who  further  generalized  (1)  by  including  an  additive  strength 
term.  Attempts  to  use  (l),  or  some  modification  of  it,  to  correlate  experimental 
values  of  U  and  have  been  extremely  successful. 

By  introducing  plausible  assumptions  concerning  the  behavior  of  the  jet  in 
flight,  Pugh"*"  was  able  to  integrate  the  penetration  velocity  and  account  for 
qualitative  features  of  the  penetration- standoff  curve.  This  method  requires 
defining  am  appropriate  penetration  velocity  U  for  the  penetration  produced  by 
the  jet  after  it  has  broken  into  a  series  of  discrete  particles.  Because  the 
correct  definition  of  U  is  not  immediately  obvious,  this  report  describes  a 
different  approach  to  the  problem  in  which  the  penetration  process  is 
represented  by  different  mathematical  models  depending  on  whether  the  jet  is 
continuous  or  broken-up. 

AH  APPROXIMATE  SOLUTION  FOB  THE  PENETRATION  VARIABLES 

To  calculate  the  total  depth  of  penetration  in  a  massive  target,  it  is 
necessary  to  add  the  penetrations  produced  by  the  various  jet  elements.  For 
a  continuous  jet,  correct  values  are  obtained  by  integrating  the  penetration 
equation.  After  the  jet  fractures  into  a  series  of  discrete  particles,  the 
physical  process  is  better  regarded  as  an  addition  of  discrete  elements  of 
penetration  corresponding  to  discrete  particles  of  jet,  and  a  more  appropriate 
mathematical  description  is  a  summation  of  finite  increments.  Accordingly, 
the  total  penetration  P^  is  given  by 


where  t  is  the  time  at  which  the  jet  reaches  the  target;  t  ,  the  time  at  which 
o  th 

the  jet  breaks  up;  and  AP^,  the  penetration  produced  by  the  I  particle. 

For  a  standoff  of  two  or  more  cone  diameters,  it  is  well  known  that  the 
velocity  of  the  forward  tip  of  the  jet  emerging  from  a  target  of  thickness 
P,  (0<  P<  P^),  can  be  calculated  (to  within  a  few  per  cent)  from  penetration- 
time  measurements  by  assuming  that  all  jet  elements  are  formed  simultaneously 
at  a  virtual  origin.  Therefore  an  approximate  solution  for  the  penetration 
variables  can  be  obtained  by  assuming  that,  at  t  =  0,  the  jet  is  formed  at  a 


virtual  origin  located  a  distance  z  in  front  of  the  target  surface.  As 
shown  in  Figure  2(a),  the  penetration  depth  and  impinging  jet  velocity  are 
related  by 

P(t)  =  tVj(t)  -  zq  ,  (3) 

where  V.(t)  is  the  velocity  of  the  jet  at  the  stagnation  point.  It  is  further 

J 

assumed  that  the  effects  of  strength  are  negligible;  hence,  for  a  continuous 
jet  (x.  =  l),  Equation  (l)  reduces  to 


U  =  Vj/(1  4  y), 


(V 


/  /  \l/2 

where  y  =  (p^/p^ j  '  .  When  (3)  is  differentiated  with  respect  to  time,  one 
obtains 


U  =  dP/dt  =  V.  +  tdV  /dt. 


(5) 


After  the  penetration  velocity  is  eliminated  between  (4)  and  (5),  the  resulting 


equation  in  V.  and  t  can  be  integrated  to  give 


Vj  =  v“(t0/t) 


,?/(!  +  y) 


where  V  A  is  the  velocity  of  the  unimpeded  foremost  tip  of  the  jet,  and  tn  is  the 
time  at  which  it  first  encounters  the  target,  relative  to  t  =  0  at  the  virtual 
origin.  Equation  (6),  together  with  (3),  can  be  used  to  compute  penetration- 
times  and  emerging- jet-velocities  once  the  constants  V°  and  tQ  have  been 
determined  experimentally,  where  Equation  (3)  is  now  given  by, 


P(t)  =V°t  (t  /t)  7^1  +  ^ 
J  O 


(7) 


After  a  time  t^  the  penetration  process  changes  from  one  involving  a 
continuous  jet  to  one  involving  a  jet  of  discrete  particles.  An  experimental 
determination  of  tn  can  be  obtained  from  the  point  at  which  the  penetration - 
time  curve  departs  from  that  predicted  for  a  continuous  jet.  If  it  is  assumed 
that  break-up  occurs  simultaneously  throughout  the  remaining  jet,  then,  as 
shown  in  Figure  2(b),  the  following  relation  must  hold: 


N 

& 

i=l 


z  +  P(t, ) 
o  v  1' 


t±vfn 


(8) 


TEC: 


01  ryi# 


nr 


where  V^n  is  "the  minimum  jet  velocity  capable  of  penetrating  the  target 


J 

material  and 


w 


Ai  is  the  length  of  the  effective  jet  that  has  not  yet 

£l  1 


penetrated  the  target.  To  calculate  the  total,  penetration,  it  is  assumed  that 
each  particle  °f  &  broken-up  jet  behaves  like  a  short  steady- state  jet"*"; 
i.e.,  AP^  =  Ai^/7.  According  to  (2)  and  (8)  and  the  definition  of  AI\,  the 
total  penetration  will  be  given  by 

PT  =  P(tx)  +  jj5o  +  p(tx)  -  t1v^inJ  //,  (9) 

which  is  independent  of  the  total  number  of  jet  particles  involved. 

Flash  radiographs  indicate  that  the  broken  jet  can  be  regarded  as  a  series 
of  discrete  particles  of  approximately  equal  length'*.  If  it  Is  assumed  that 


Ai.  =  Ai_  = 
1  d 


.  Ai  =  Ai.  the  number  of  effective  jet  particles  will  be  given  by 


N 


=  [z0  +  p^)  -  ty^n  /Ai. 


(10) 


Flash  radiographs  also  show  that  the  increment  in  velocity  between  adjacent 
particles  is  approximately  a  constant-7.  Hence,  the  increment  in  velocity 
between  adjacent  particles  is  given  by 


av  =  -  Tv  ft.  ^  -  v^11 

L'j''1'  J 


/fw  _  1  1 


fin 

x-1-— / 


or  1  +"h*=»  Trml  aat  4"ir  rcP  ■? 

uut  v  -w'-xuj  ui  out  J_ 


.  th 


T~»C1  -V>4-  ^  rtl  , 

^ic*-L . 


•?  r i 

-lo  51  V  Cil  uj 


(bb  -  VA>  +  h  -  ^V 


fl2l 


An  inspection  of  Figure  2(c)  shows  that  the  time  between  the  arrival  of 

, ,  .  th  .  , ,  / .  .  %  st  ... 

xne  1  ana  xne  (1  t  _l)  parxicie  xs  given  Dy 


Ati,i+i  =  +  “  -  AVj(ti  -  h>]  /<vjW  (15) 


Of  course,  the  arrival  time  of  the  particle  is  given  by 

i 


: .  =  tn  t  \  At  .  .  , 

i  1/  J-1,0 


(It) 


1-0 


where  i  >  2. 
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The  penetration  depth  at  the  instant  the  i  particle  reaches  the  bottom 
of  the  hole  is  given  by 

P(tj_)  =  P(tx)  +  (i  -  1)AP.  (15) 

The  penetration  by  the  discontinuous  portion  of  the  jet  will  be  completely 
specified  once  the  constants  Ai  and  V?111  have  been  determined.  The  number  of 
particles  is  determined  from  (10),  and  their  respective  velocities  are 
determined  from  (ll)  and  (12).  The  time  at  which  each  particle  starts  to 
penetrate  the  target  can  be  determined  from  (13)  and  (14),  and  the  corresponding 
penetration  depths  can  be  obtained  from  (15) • 

The  accuracy  of  the  theory  for  correlating  the  penetration  variables  is 

best  illustrated  by  analyzing  data  obtained  with  the  105  nun  unconfined  test 

charge  used  at  the  Ballistic  Research  Laboratories.  The  time  at  which  fracture 

of  the  jet  occurs  was  determined  by  integrating  the  hydrodynamic  penetration 

equation  for  a  continuous  jet  and  comparing  it  with  the  experimental  penetration- 
* 

time  data  .  The  point  at  which  the  two  diverge;  i.e. ,  the  time  when  the 

continuous  jet  approximation  no  longer  adequately  describes  the  penetration 

process,  is  taken  to  be  the  time  at  which  simultaneous  fracture  occurs  in  the 

remainder  of  the  effective  jet.  As  shown  in  Figure  3>  the  divergence  of  the 

experimental  data  from  the  continuous  jet  approximation  occurs  at  120  p-sec. 

The  value  of  Ai  was  estimated  from  flash  radiographs  to  be  1.68  cm.  The  values 

of  V™111  and  V°  were  measured  directly  and  found  to  be  0.210  cm/ii-sec  and 
<J  J 

0.701  cm/p-sec,  respectively. 

The  penetration  theory  for  discrete  particles  can  be  used  to  extend  the 
penetration-time  curve  beyond  120  p-sec.  As  shown  in  Figure  3,  results  of  the 
calculation  are  in  good  agreement  with  experimental  observations.  In  addition 
to  the  penetration-time  data,  the  velocity  of  jets  emerging  from  targets  of 
various  thicknesses  gives  an  independent  measure  of  the  jet  velocity  as  a 
function  of  the  penetration  depth.  As  illustrated  in  Figure  4,  the  curve 
computed  from  the  penetration  theory  adequately  represents  the  emerging  jet 
velocity  measurements. 


See  reference  4,  for  a  general  description  of  the  experimental  techniques 
used  in  obtaining  penetration -time  data. 
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'THE  USE  OF  CONTINUOUS  VARIABLES  TO  REPRESENT 
THE  PENETRATION  PRODUCED  BY  A  PARTICLE  JET 

It  has  already  been  noted  that  the  total  penetration  produced  by  the 
broken-up  portion  of  the  jet  depends  on  the  total  length  of  all  the  particles 
involved,  but  not  on  the  length  of  an  individual  particle.  It  is  therefore 
permissible  to  represent  the  penetration  variables  by  continuous  functions 
obtained  from  appropriate  mathematical  limits  taken  as  the  particle  length  AZ 
approaches  zero  and  the  number  of  particles  N  approaches  infinity.  Of  course- 
the  limits  must  be  taken  in  such  a  way  that 

L  =  Lim  (I m)  =  +  P(tJ  -  Vv™-n  (lb) 

Ai-o  0  1  1  J 

N  — » * 

is  a  constant. 


An  expression  for  the  velocity  of  the  impinging  jet  is  obtained  after 
dividing  (ll)  by  AZ, 


AV./AZ  =  -  [VjCy  -  vfn]  /(mi  -  az), 

taking  the  limit  as  AZ  tends  to  zero. 

Urn  (AV./Al)  .  dV  .jit,  >  -  fv.(t.  )  -  V-j  /L 
Ai-,0  J  J  UO'  I'  j  J  ■ 

and  integrating, 

^Vj(s)  s 

dV.  =  -  (i/L)  |,.x„i 

A) 


V  .  (t  )  -  v1^ 


n 

nin  [ 

’  Vo 


a  e 

U.  x/  • 


(17) 


( iA\ 
\^  ) 


(19) 


The  final  expression  for  the  jet  velocity  is 


V.(s)  =  V  ,  (t. )  -  fv  .(t.)  - 

j  j  i  L  J  1  j  . 

f 


s/L, 


where  the  parameter  s 


d Z  is  the  location  of  an  infinitesimal 


element  of  jet  relative  to  the  stagnation  point  at  time  t1* 


(20) 


12 


An  expression  for  the  penetration  is  readily  obtained  by  taking  the  limit 
of  AP^  as  Ai  tends  to  zero, 

Lim  {AP/A  l)  =  dP/di  =  l/y  (21) 

Ai  -*0 


aud  integrating. 


P(s) 

far- 


p(h) 


,,r, 

i 

Jo 


ft  l..\  I  J  * 

k-L/7;  I  n*. 


^  Cm(—  f 


The  resulting  expression  for  the  penetration  is 
P(s)  =  P(t1)  +  s/7. 


(25) 


In  much  the  same  way,  an  expression  for  the  time  is  obtained  after 
dividing  (13)  by  AJ. 

(Vj)i+l(Ati,i+l/Ai)  =  (**/&*)  +  1  -  (A  Vj/Ai)  (t.  -  t1)  ,  (2k) 

taking  the  limit  as  AZ  tends  to  zero. 


V.dt/di  =  (1/7)  +  1  -  (t.  -  tJdV./di  , 


(25) 


and  integrating 


r 


v 


d(V.t) 


J 

rr  t  x.  \  4- 

vjvvti 


(1  +  r) 

7  7. 


di  +  t 


.  r 


v,(s) 


dV  .. 

j-  j 

ir  /+  \ 
Vj'“l' 


(26) 


The  final  expression  for  the  time  is 


t(£ 


~  V  / 


i)M.(s). 


(27) 


Equations  (20),  (23),  and  (27)  are  parametric  equations  from  which  the 
penetration  variables  can  be  calculated.  The  parameter  is  permitted  to  vary 
continuously  over  the  range  0  <  s  <  L,  and  the  penetration  variables  calculated 
from  the  equations  are  continuous  functions.  A  graph  illustrating  the 
functional  relation  between  the  penetration  and  the  time  is  presented  in  Figure  5* 
For  comparison  purposes,  results  computed  from  the  discrete  particle  theory 
have  been  superimposed.  The  time  and  penetration  depth  corresponding  to  each 


particle  impact  on  the  target  were  computed  from  (l4)  and  (15)  respectively. 

To  emphasize  the  discrete  nature  of  the  process,  the  time  required  for  each 
particle  to  complete  its  penetration  was  computed  from 

At.  =  AP./U.  =  (1  +  7)  P^V  .).  •  (28) 

According  to  the  idealized  model  of  the  particle  jet,  there  will  be  intervals 

,  \  * 
of  time  (At^  -  At^)  during  which  no  penetration  takes  place  .  The  length 

of  this  interval  increases  with  the  depth  of  penetration  because  adjacent  particles 

have  more  time  to  spread  due  to  differences  in  their  velocities.  As  shown  in 

Figure  5,  the  continuous  parameter  representation  of  the  penetration-time  curve 

follows  the  discrete  particle  theory  very  closely.  Because  the  equations  are 

simpler,  future  calculations  of  penetration  variables  will  use  the  continuous 

parameter  representation. 


For  a  continuous  jet  Bernoulli's  equation  was  used  to  establish  a  relation 
between  the  two  velocities  V.  and  IT*".  A  similar  relation  can  be  obtained  for  a 

d 

jet  of  discrete  particles  by  differentiating  (20),  (23)  and  (27)  with  respect  to 


t ,  obtaining 


dV 


j/dt  -  -[v.^)  -  v^nJ  (ds/dt)/L  , 
U  =  dP/dt  =  (1/7)  ds/dt  , 


(29) 

(30) 


(t  -  tx)dV  /dt  +  Vj  =  £(7  +  D/t]  ds/dt  (3l) 

respectively.  When  (30 )  is  used  to  eliminate  ds/dt  from  (29)  and  (31),  one 
obtains , 

dVj/dt  =  -  [vj(tl}  -  v^]  7u/l  ,  (32) 

(t  -  t,)dV./dt  +  V.  =  (7  +  l)U.  (33) 


The  discontinuous  nature  of  the  penetration  variables  was  neglected  in 
Figure  3  and  Figure  4,  a  smooth  curve  being  drawn  through  points  representing 
the  variable  at  the  instant  each  particle  impacts  with  the  target. 
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Elimination  of  dV./dt  between  (32)  and  (33)  gives 
J 

-  [vjCt^)  -  vfn]  7U(t  -  tx)/L  +  Vj  =  (7  +  1)U. 


(34) 


To  eliminate  (t  -  t1 )  from  the  preceding  equations,  s  is  eliminated  between 
(20)  and  (27)  to  give 


fx. 


tx)  =  (7  +  1)L 


.(t-i)  - 


1 


/».t! 


■j'-v  -  uj  "’j  l’j'v 

Finally,  substitution  of  (35)  into  (34)  gives 


Trmin“ 

V  . 

J 


f^"\ 
\ss  ) 


U  =  vy<7  +  l)V.(tn). 


(56) 


For  a  jet  of  discrete  particles,  the  penetration  velocity  is  proportional  to 

the  square  of  the  jet  velocity,  and  for  this  reason,  the  properties  of  a 

discontinuous  jet  cannot  be  described  by  replacing  p  in  Equation  (l)  with  an 

J 

effective  density.  It  is  interesting  to  note  that  the  problem  of  defining  an 
appropriate  penetration  velocity  for  a  discontinuous  jet  has  been  solved,  and 
the  penetration  can  be  obtained  by  integrating  the  penetration  velocity;  i.e., 

rt  Al  A 

b(t)  =  J  Udt  =  j  U^at  +  j  Ugdt,  (t  >  t1).  (37) 

U  L/  Vj 

00  1 

Where  U,  is  given  by  (4)  and  UU  by  (3 6). 

JL  w  ”  x  ^  ‘  ' 


SCALING  RELATIONS 

It  has  been  shown  experimentally  that  the  penetration-time  curves  for  a 

homologous  series  of  scaled  charges  reduce  to  a  single  curve  when  both  the 

5 

penetration  and  time  are  divided  by  the  charge  diameter  .  The  accuracy  with 
which  a  single  curve  represents  the  penetration- time  data  is  illustrated  in 
Figure  6.  For  the  data  presented  in  Figure  6,  the  location  of  the  virtual 
origin  with  respect  to  the  top  of  the  target  can  be  expressed  as  z Q  =  t 
where  D  is  the  charge  diameter  and  £  ,  a  constant. 

Equations  (3)  and  (6)  can  be  used  to  show  that  the  penetration,  measured 
in  charge  diameters,  is  given  by 

■p  /r>  _  rr\  —  fry  -  !"  f 

r/ii  •-  -  iy^  -  i,o  ,  w  w 
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where 


V .  =  V  ,°  (t  +  y) 

J  j  o'  ' 


(59) 


and  t  =  t/D.  (40) 

After  the  jet  breaks  up  the  penetration  can  be  expressed  as  a  function  of  the 
dimensionless  parameter  a  =  s/d.  The  penetration,  measured  in  charge  diameters, 
is  given  by 

p/D  =  jt(cf)  =  rt(i1)  +  a/y  (4-1) 

and  the  corresponding  times  are  given  by 


t/D  =  t(c)  =  x1  +  0(7  +  l)/rV,  , 


(42) 


where  V  .  is  now  given  by 


“1 

r 

r..  f  ,  „  min 

Fax.;  -  V  . 

J  1  0 

L_  — 1 

«/  U.  ♦  -  TiviMn 

1—  ”  — 

r\,*\ 
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The  preceding  equations,  (38)  to  (43),  describe  the  penetration  of  a  shaped- 

charge  jet  in  terms  of  size-  independent  variables,  x  =  t/D  and  n  =  p/D,  and 

are  obtained  from  Equations  (l6),  (20),  (23)  and  (27).  For  a  homologous 

series  of  scaled  charges  there  exists  a  unique  relation  between  «  and  x,  which 

is  independent  of  the  charge  size.  The  penetration  velocity  U  and  the  jet 

velocity  V  .  can  also  be  represented  as  functions  of  either  n  or  x,  the  functions 
J 

being  independent  of  charge  size. 


GENERALIZATION  OF  THE  THEORY  FOE  SHORT  STANDOFF 


The  present  theory  will  be  inaccurate  for  short  standoffs  where  it  is  no 
longer  permissible  to  neglect  small  differences  between  the  times  and  positions 


at  which  different  jet  elements  are  formed.  To  remove  the  inaccuracies  at 
short  standoffs,  a  complete  description  of  the  liner-collapse  and  jet- formation 
process  is  required.  At  the  present  time  the  necessary  data  are  available  for 


only  one  charge,  the  105  ran  unconfined  test  charge  used  at  the  Ballistic 

6 

Research  Laboratories  .  For  this  charge,  the  collapse  velocity  in  cm/p-sec 
imparted  to  an  element  of  the  liner  is  given  by 
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V 

o 


=-0.197898  +  0.288477x-  0.0775152x2  +0.009426l5x5 
-  0.00044 5364x\  for  4.0  cm  <  x  <  8.0  cm. 


m 


where  x  is  the  axial  position  of  the  liner  element  measured  relative  to  the 
apex  of  the  conical  liner. 


Equations  describing  a  point  initiated  charge  of  length  L  containing  a 
conical  cavity  of  height  H  are  readily  derived  using  the  quantities  illustrated 
in  Figure  7*  The  zero  reference  for  time,  t  =  0,  is  taken  to  he  the  instant 
the  explosive  is  initiated.  The  time  required  for  the  detonation  wave  to  reach 
a  liner  element  located  at  the  position  x  is  given  by 

T(x)  =  |"(L  -  H  +  x)2  +  (x  tan  a)2"]  1/^2/uri  ,  (45) 


where  is  the  detonation  velocity  of  the  explosive.  The  time  at  which  the 
liner  element  reaches  the  axis  is  obtained  from 

t+  (x)  =  T(x)  +  x  tan  a/Vocos  (a  +  5)  (46) 

and  the  point  at  which  the  liner  element  reaches  the  axis  is  given  by 


z+(x)  =  x  +  (t+  -  T)VQsin  (a  +  5) . 

Values  for  the  angle  5  can  be  computed  from 

5  =  sin  ^(Vocos e  /2Up), 

where  tan(90°  -  a  +  e)  =  (L  -  H  +  x)/(x  tan  a). 

The  depth  of  penetration  at  a  time  t  will  be  given  by 

\  /j_  J--N.T  In  -r  _+\ 

r\T.)  =  (.u  -  -g  +  a  -  z  ;  , 


(47) 

(48) 

(49) 


/ 1-^\ 


where  S  is  the  standoff  between  the  base  of  the  charge  and  the  top  of  the  target. 

The  jet  velocity,  V.(t),  appearing  in  the  penetration  equation  is  the  value  for 

the  jet  element  at  the  stagnation  point  between  the  jet  and  target.  The  jet 

velocity  may  also  be  regarded  as  a  function  of  x,  the  original  position  of  the 

element  within  the  liner.  The  function  V.(x)  can  be  computed  from  the  non-steady 

2  J 

theory  of  jet  formation  .  The  penetration  velocity  obtained  by  differentiating 
(50)  is 

U  =  dP/dt  =  Vj  +  -  Q(Vj)J  (dVj/dt)  ,  (51) 
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where 


, .  +  / 


+  / , 


Q (V j j  =  (f  +  Vj  df/dV.  -  dz/dV.j.  (52) 

The  differential  equation  obtained  by  eliminating  U  between  (51)  and  (4)  is 
dV  ./dt  =  -yV  ./ (l  +  7)  ft  -  Q(V,)1  .  (55) 

J  j  I —  J  — 1 

The  virtual  origin  approximation  is  obtained  by  letting  Q(V.)  be  zero. 

J  t 

Equation  (53)  is  a  first  order  linear  differential  equation  of  the  form  dV./dt  = 

f (V  .,t)  and  can  be  solved  by  well  known  numerical  procedures. 

J 

The  numerical  solution  to  the  differential  equation  (53)  requires  a  precise 
knowledge  of  the  cone  collapse  parameters.  The  data  presently  available  on  the 
105mm  unconfined  test  charge,  while  adequately  describing  the  jet  formation,  do 
not  provide  sufficiently  precise  initial  conditions  necessary  to  solve  the  short 
standoff  equation.  In  lieu  of  using  the  generalization,  the  virtual  origin 
approximation  was  shown  to  be  satisfactory  for  standoffs  down  to  half  of  a  cone 
diameter.  Hence  the  generalization  reduces  to  one  of  academic  interest,  and  will 
be  resolved  when  more  precise  data  becomes  available. 

EFFECTS  OF  COMPRESSIBILITY 

As  pointed  out  in  the  introduction,  the  penetration  of  a  target  by  a 


shapect  charge  jen  is  aescrioea  oy  incompressiu-Le  iiuiu  now.  nuwever, 
compressibilities  of  the  target  and  jet  are  certainly  not  negligible  at  the 
pressures  produced  during  the  penetration,  and  it  is  important  to  estimate 
the  error  introduced  by  the  incompressible  approximation. 

Consider  the  steady  state  penetration  of  a  compressible  target  by  a 
compressible  jet  in  which  the  jet  and  target  materials  approach  the  stagnatior 
point  with  subsonic  velocities.  For  adiabatic  compressible  flow,  Bernoulli's 
equation  requires  that 


I 


dp  -  (l/2)p°U2  =  0, 


(54) 


o. 

"t 


^  dp  -(1/2)  P°(V.  -  U)2 
Pj 


(55) 
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along  the  axial  streamline,  where  pg  is  the  stagnation  pressure ;  and  p°,  the 
respective  densities  at  normal  pressure.  The  effect  of  compressibility  on  the 


penetration  velocity  is  readily  determined  by  setting  the  left-hand  side  of  (54) 
equal  to  the  left-hand  side  of  (55)  and  rearranging  terms.  Thus  one  obtains 


(1/2) p°  (V j 


U)2  -  (1/2)  p°\f= 


(p>j)  -  (Pt°/Pt) 


dp. 


(56) 


If  the  material  is  incompressible  the  integral  in  (56)  is  zero  and  the  resulting 
expression  is  equivalent  to  (l)  with  X  set  equal  to  one.  It  is  apparent  that  the 
integral  is  zero  if  the  jet  and  target  are  the  same  material.  On  the  other  hand, 
the  largest  error  will  result  from  treating  as  incompressible  two  materials  that 
have  greatly  different  compressibilities. 


It  is  possible  to  compute  the  error  for  any  combination  of  materials  for 

which  suitable  equation-of -state  data  are  available.  Usually,  the  high  pressure 

dynamic  measurements  provide  an  empirical  relation  between  the  pressure  and  density 

along  the  shock  Hugoniot,  When  the  Gruneisen  ratio  is  also  known,  the  pressure 

and  density  can  be  calculated  for  any  adiabat  that  intersects  the  Hugoniot+. 

The  penetration  of  a  titanium  target  by  a  copper  jet  was  chosen  as  an 

illustration  because  titanium  is  considerably  more  compressible  than  copper. 

Calculations  were  made  using  the  equation  of  state  data  reported  by  McQueen  and 
v 

Marsh1 .  Assuming  a  stagnation  pressure  of  500  kilobars ,  a  penetration  velocity 

of  0.412  em/p-sec  is  computed  from  (54),  and  a  jet  velocity  of  O.726  cm/p-sec  is 

obtained  from  (55)-  For  a  penetration  velocity  of  0.412  cm/p-sec,  the 

incompressible  approximation  predicts  a  jet  velocity  of  0.720  cm/p-sec,  which 

is  within  1%  of  the  more  exact  calculations.  The  relative  magnitudes  of  the 

quantities  appearing  in  (56)  are  illustrated  graphically  in  Figure  8.  The 

between  (l/2  )p°(V.  -  U)2  and  (l/2  )p°TJ2,  which  is  represented  by  the 
J  J  u 


For  a  general  discussion  concerning  the  behavior  of  solids  at  high  dynamic 
pressures,  the  reader  is  referred  to  the  section  entitled  "Compression  of 
Solids  by  Strong  Shock  Waves"  contributed  by  M.  H.  Rice,  R.  G.  McQueen,  and 
J.  M.  Walsh  to  Solid  State  Physics,  Vol.  6,  Edited  by  Fredrick  Seitz  and 
David  Turnbull,  Academic  Press,  Inc.  Publishers  (1958). 
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shaded  area,  is  a  small  fraction  of  ( l/ 2 )  p°u  ,  which  is  represented  by  the 
cross  hatched  area.  It  is  noted,  however,  that  the  pressure  computed  from  the 
incompressible  approximation,  (l/2)  p°U^,  is  422  kilobars,  or  1 6$  less  than  the 
actual  pressure.  Thus,  as  already  demonstrated  by  direct  experimental 
observations,  the  relationship  between  the  jet  and  penetration  velocities  obtained 
from  the  incompressible  approximation  is  sufficiently  accurate  for  all  practical 
problems  involving  the  penetration  of  shaped-charge  jets . 
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DURING  JET  FORMATION 


DURING  PENETRATION 


Fig.  1  Line  drawings  illustrating  the  formation  of  a  shaped- 
charge  jet  and  the  subsequent  penetration  of  a  target 
by  the  jet.  Note  the  hydrodynamic  nature  of  the 
penetration  process. 


VIRTIIAl  ORIGIN 
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virtual  origin  approximation 


penetration  variables 
and  (c)  the  time  interval 
i  the  (i  +  l)st  particles. 


TIME  (  M  SEC  ) 

Fig.  3  Penetration-time  curve  for  the  105mm  unconfined  test  charge  fired  at  an 
8  1/4  in.  standoff.  Note  that  the  continuous  jet  approximation  is  valid 
lor  the  first  120p-see.  Beyond  this  point,  the  penetration  time  must 
be  calculated  using  a  method  applicable  to  a  jet  of  discrete  particles. 


0  4  8  12  16  20 

PENETRATION  (IN) 

Fig.  4  Emerging  jet  velocity  versus  the  penetration  depth  for 
the  105mm  unconfined  test  charge  fired  at  an  8  1/4  in. 
standoff.  The  theoretical  curve  was  computed  from  the 
penetration  theory  developed  in  this  report. 


25 


120  140  160  180  200  220  240  260  280  300 

TIME  (/A  SEC  ) 

Fig.  5  Curves  representing  the  penetration  produced  by  the 
broken-up  portion  of  jets  obtained  from  105mm 
unconfined  test  charges  fired  at  an  8  l/b  in.  standoff. 


PENETRATION  IN  CONE  DIAMETERS 


T  IN  ^iSEC./CM 


Fig.  6  Graph  showing  that  the  penetration-time  data  for  a  homologous  series  of  scaled  charges  can. 
be  represented  by  a  single  curve  when  both  the  penetration  ang.  time  are  divided  by  a 
characteristic  linear  dimension,  such  as  the  charge  diameter.-5 
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Fig.  8  Graphical  illustration  of  the  relative  magnitudes  of  the 

quantities  resulting  from  the  compressible  and  incompressible 
treatment  of  Ti  penetrating  Cu. 
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In  earlier  attempts  to  describe  the  penetration  of  shaped-charge  Jets,  the 
total  depth  of  penetration  was  computed  from  J  out.  Evaluation  of  the  integral 
has  not  been  possible  because  the  penetration  velocity  U  for  the  discontinuous 
portion  of  the  jet  has  not  been  properly  defined.  It  is  nov  postulated  that  the 
[total  penetration  can  he  represented  by  f  Udt  +  ^~’AP,  where  the  integral  is  taken 
over  the  penetration  produced  uy  the  CuutlnUOUB  portion  Ox  the  jet  Slid  the 
5 1  lrrnflt. Inn  j_R  ~hflkp>n  over  the  penet-ratinn  produced  by  the  discrete  parti cies.  It  is 
i shown  that  numerical  values  of  the  penetration  variables  for  the  105™n  unconfined 
i charge  are  compatible  with  the  modified  penetration  equation.  Scaling  relations 
are  introduced,  and  the  penetration  equations  are  given  In  size  independent 
'variables  Although  it  was  assuz^ed  that  the  Incompressible  approximation  provides 
an  adequate  relationship  between  the  Jet  and  penetration  velocities  for  a 
continuous  Jet,  it  is  later  shown  that  the  finite  compressibilities  of  the  Jet  and 
'target  materials  do  not  seriously  alter  this  relation,  xhe  incompressible 
approximation  does  underestimate  the  pressure  by  an  appreciable  amount. 
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,  in  earner  attempts  to  aescriDe  tne  penetration  or  snapea-cnarge  jets,  the 
tot  el  depth  of  penetration  was  computed  from  J* Udt.  Evaluation  of  the  integral 
tias  not  been  possible  because  the  penetration  velocity  U  for  the  discontinuous 
portion  of  the  Jet  has  not  been  properly  defined.  It  ia  now  postulated  that  the 
total  penetration  can  be  represented  by  j  Udt  4  V  AP,  where  the  integral  is  taken 
I  over  the  penetration  produced  by  the  continuous  portion  of  the  Jet  sad  the 
'summation  is  taken  over  the  penetration  produced  by  the  discrete  particles.  It  Is 
{Shown  that  numerical  values  of  the  penetration  variables  for  the  103m  unconfined 
i  charge  are  compatible  with  the  modified  penetration  equation.  Scaling  relations 
are  introduced,  and  the  penetration  equations  are  given  in  size  independent 
! variables.  Although  it  was  assumed  that  the  Incompressible  approximation  provides 
|  an  adequate  relationship  between  the  jet  and  penetration  velocities  for  a 
{Continuous  jet,  it  is  later  shown  that  the  finite  compressibilities  of  the  jet  and 
target  materials  do  not  seriously  alter  this  relation.  The  incompressible 
{approximation  does  underestimate  the  pressure  by  an  appreciable  amount. 
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In  earlier  attempts  to  describe  the  penetration  of  shaped-charge  jets,  the 
total  depth  of  penetration  was  computed  from  J'  Udt.  Evaluation  of  the  integral 
has  not  been  possible  because  the  penetration  velocity  U  for  the  discontinuous 
j  portion  of  the  jet  has  not  been  properly  defined.  It  is  now  postulated  that  the 
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approximation  does  underestimate  the  pressure  by  an  appreciable  amount. 
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In  earlier  attempts  to  describe  the  penetration  of  shaped-charge  jets,  the 
total  depth  of  penetration  was  computed  from  j  Udt .  Eval'intion  of  the  Integral 
has  not  been  possible  because  the  penetration  velocity  U  for  the  discontinuous 
portion  of  the  Jet  has  not  been  properly  defined.  It  is  now  postulated  that  the  . 
total  penetration  can  be  represented  by  J~  Udt  +  J^AP,  where  the  Integral  is  taken 
over  the  penetration  produced  hy  the  continuous  portion  of  the  jet  the 
summation  is  taken  over  the  penetration  produced  by  the  discrete  particles.  It  la 
shown  that  numerical  values  of  the  penetration  variables  for  the  105nmi  unconfined 
charge  are  compatible  with  the  modified  penetration  equation.  Scaling  relations 
are  introduced,  and  the  penetration  equations  are  given  in  size  Independent 
variables.  Although  it  was  assumed  that  the  incompressible  approximation  provides 
an  adequate  relationship  between  the  jet  and  penetration  velocities  for  a 
continuous  Jet,  it  is  later  shown  that  the  finite  compressibilities  of  the  Jet  and 
target  materials  do  not  seriously  alter  this  relation.  The  incompressible 
1  approximation  does  underestimate  the  pressure  by  an  appreciable  amount. 


